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ABSTRACT 

We study the scattering of an incident wave by a flexible panel. The panel vibration 
is governed by the nonlinear plate equations while the loading on the panel, which is the 
pressure difference across the panel, depends on the reflected and transmitted waves. Two 
models are used to calculate this structural-acoustic, interaction problem. One solves the 
three dimensional nonlinear Euler equations for the flow-field coupled with the plate equa- 
tions (the fully coupled model). The second uses the linear wave equation for the acoustic 
field and expresses the load as a double integral involving the panel oscillation (the decoupled 
model). The panel oscillation governed by a system of integro-differential equations is solved 
numerically and the acoustic field is then defined by an explicit formula. Numerical results 
are obtained using the two models for linear and nonlinear panel vibrations. The predictions 
given by these two models are in good agreement but the computational time needed for the 
“fully coupled model” is 60 times longer than that for “the decoupled model”. 


1 This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract Nos. NAS 1-19480 and NAS 1-18605 while the third author was in residence at the Institute for Com- 
puter Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 
23681-0001. 
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1. Introduction 


We consider the scattering of an incident acoustic wave by a planar interface, the (;r, 
z) plane. The interface is a rigid surface except on the flexible thin panel, T>, with panel 
thickness h much smaller than the size L of the panel, or the size of Z>, 

h/L< 1. (1.1) 

This is a basic structural-acoustic interaction problem simulating the transmission of an 
external acoustic source, say the engine noise, through the airframe into its interior. It 
is an interaction problem because the incident wave excites panel oscillation which in 
turn induces transmitted and diffracted waves. Meanwhile, these waves contribute to the 
loading on the panel. Although in an engineering problem, the typical amplitude of the 
transverse displacement, ?/* , has to be much smaller than the size of the panel, i. e., 

V*/L<. 1 , ( 1 . 2 ) 

the problem can be nonlinear in three different aspects: 

(1) The panel oscillation can be nonlinear when the transverse displacement of the panel 
is not much smaller than the panel thickness h, that is 

77* ss h , ( 1 . 3 ) 

Nonlinear plate equations are needed usually when the incident wave is in resonance 
with the panel oscillation. 

( 2 ) In the near-field, of length scale L. the flow-field is linear and obeys the simple wave 
equation but in the far-field, of length scale much larger than L, the flow-field may 
become nonlinear when the second order terms are needed to account for the gradual 
steepening of compression waves. 

( 3 ) I11 case the flow-field is nonuniform and/or the initial pressure variation is no longer 
much smaller than the ambient pressure then we have nonlinear flow in the near-field. 
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We have to solve the Euler equations with appropriate initial data to simulate an 

incident wave and far-field conditions. 

Regardless of whether the flow-field is linear or nonlinear, the plate should have sufficient 
rigidity so that the transverse displacement of the plate remains small and hence the 
governing equations for the plate derived under assumption (1.1) remain valid. 

For the problem of sound transmission through a moving airframe, the external acous- 
tic field is governed by the convective wave equation while that in the interior by the simple 
wave equation. Also the pressure and density in the exterior can be much lower than those 
in the interior. In the simple model simulating the experimental data with no-flow on both 
sides of the panel, both acoustic fields are governed by the simple wave equation with the 
same ambient pressure and density. The validity of the linear theory for the flow-field is 
confirmed by the investigations of Frendi et al. [1,2] and Maestrello et al. [3]. They showed 
that even when the panel vibration is nonlinear, the acoustic field can remain linear. 

The pressure difference across the panel induces the panel oscillation which in turn 
excites the transmitted wave and an additional reflected wave. Thus the solution of the 
panel oscillation is coupled with the solutions of the scattered and transmitted waves which 
satisfy the far-field radiation conditions for outgoing plane waves (see for example [4]). In 
general numerical solution of the scattering of an incident wave by an elastic scatterer, 
an interface, or a panel, requires the introduction of a finite computational domain for 
both acoustic fields. Higher order radiation conditions were derived [5-6] so that they can 
be imposed on the boundary of the finite computational domain to give more accurate 
approximation to the solution in the unbounded domain. Since the size of the computa- 
tional domain has to be much larger than the size of the panel or a scatterer, numerical 
solution of this three-dimensional unsteady problem is very tedious especially when the 
computation has to be continued for a long time relative to the period of oscillation of the 
panel. Furthermore, the accuracy of the solution depends not only on that of the numeri- 
cal solution of the differential equations but also on the approximate boundary conditions. 
Refinements of the grid size and time step have to go with an enlargement of the computa- 
tional domain. It is desirable to find exact boundary conditions so that the computational 
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domain does not have to be much larger than the scatterer and is independent of the choice 
of grid size. The exact conditions were presented in [7] for a scatterer or inhomogeneous 
medium of bounded support. With the scatterer inside the computational domain, the 
integral representation of the solution of the simple wave equation is applicable to the 
region outside and on the boundary of the computational domain. 

Here the scatterer is a panel embedded in a rigid plane. The integral representation 
can be applied on the panel and the integrand involves only the normal velocity of the 
panel [8, 9]. Thus we have a system of integro- differential equations for the panel oscillation 
decoupled from the acoustic field. This system is referred to as the decoupled model. 

For the case that the acoustic wave length is much smaller than L, the integral repre- 
sentation can be approximated by the local derivatives of the panel oscillation and hence 
the integro-differential equations are reduced to partial differential equations [8, 9]. This 
approximation is not applicable here because in the experimental data, the acoustic wave 
length is comparable to L. We need numerical solution of the integro-differential equations 
for the panel oscillation and then we can use the integral representation to evaluate the 
acoustic field and the nonlinear waves in the far-field by Whitham’s theory or matched 
asymptotics. 

Note that the solution of the decoupled panel oscillation is an unsteady two dimen- 
sional problem in the finite domain T> and is several orders of magnitude simpler than the 
solution of the fully coupled problem. To show the efficiency of the decoupled model and 
the accuracy of the solution even when the panel oscillation is nonlinear, we compare the 
solution of this decoupled system with the solution of the fully coupled problem for which 
we solve the three-dimensional nonlinear Euler equations for the flow-field coupled with 
the nonlinear equations for the panel oscillations. 

In the next section we present the mathematical formulations of these two models and 
discuss their applications in general. Section 3 describes the various numerical techniques 
used to solve the problems involved. The results and discussion are given in section 4, and 
the conclusions are in section 5. 
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2. Formulation of the Two Models 


The physical problem being studied experimentally is that of the nonlinear oscillations 
of a typical aircraft panel excited by harmonic plane waves at normal incidence. The panel 
is clamped onto a large rigid plate. Let the typical panel be represented by a rectangular 
domain V in the (x,z) plane, with length L, width W and one vertex located at (xo, ~o), 
i. e., 

V = {(x, z) \ xo < x < xo + L , zo < z < zo + W) . (2-1) 

The compliment of V in the (x,z) plane represents the rigid plate. The incident side is 
the half space y > 0 and the transmitted side is the half space y < 0. Both sides have 
the same ambient pressure p a Q and temperature Too and hence the same speed of sound C 
and density poo- We denote the pressure and velocity potential of the flow-field by p and 
<f> and use the superscripts + and — to denote the quantities on the sides y > 0 and y < 0. 
On the side y > 0, the superscripts i and r are used to denote the quantities associated 
with the incident and reflected waves, respectively. These two waves are mirror images 
with respect to the (x, z) plane. When the incident waves are plane waves advancing in 
the direction opposite to the y axis and hitting the (x, z) plane at, t 0, we hav< 

<p^ l \t, x, y, z) = f(Ct + y) and <j>^ r \t,x,y, z) = f(Ct — y) , (2-2) 

with /(£) = 0 for £ < 0. Therefore =0 when 0 > y < Ct. 

Transverse oscillation, i] ■ of the rectangular flexible panel, T), is excited by the pressure 
difference across the panel for t > 0 and the oscillation in turn induces scattered waves 6° 
in the incident side and transmitted waves 4> l . Both waves have the homogeneous initial 
data at / = 0, 

^ = 0 ’ V > ° ’ (2.3) 

4> {t) = 4>\ l) = 0 , y < 0 . 

Under linear theory, we have 

</> (r) (t,x,y,z) = 4> (0 (t,x,-y,z) ^ 

(f> W {t, x , -y, z) = -^ (s) (t, x, y, z) for y > 0 . 
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The wave fronts of <f > ^ and <f>^ are the envelopes of the sonic spheres with radius Ct and 
centered on the panel in y > 0 and y < 0, respectively. The pressure in the incident and 
transmitted sides are 


p + - Poo = p (,) + p (r) + p (a) i y > o 
p~ - Poo = p (<) , y < o 

and the pressure difference across the panel is 

Ap = p~ — p + = 2p^ - + p* r 'j , y = 0 , (x,z)€V. 


(2.5) 


(2.6) 


In this section we formulate two mathematical models to compare with and compli- 
ment the experimental studies. We present the system of equations governing the nonlinear 
flow-field and panel oscillation (the fully coupled model) in Section 2.1 and formulate the 
system of integro-differential equations for nonlinear panel oscillation including the effect 
of a linear acoustic field (the decoupled model) in Section 2.2. 

2.1 The Fully Coupled Model 

In this model the structural-acoustic interaction is analyzed by solving the three- 
dimensional nonlinear Euler equations together with the nonlinear plate equations. The 
configuration of the computational domain is shown in Fig, 1. In cartesian coordinates, x, 
y and z , the compressible, nonlinear Euler equations can be written in conservation form 
as 

Qt = F X + G y + H z , (2.7) 

where Q is the vector (p, pu , pv , pw , e), p is the density, pu , pv and pw are the x, y and 
z momenta respectively, and e is the total energy per unit volume given by 

e = ]^p(u 2 + v 2 + w 2 ) + pc v T, (2.8) 

with c v being the specific heat at constant volume. In eq. (2.7), the functions F, G and 



H are: 


F = 


/ pu \ 
pu 2 + p 

puv 

puw 

\u(e + p) } 


G — 


( pV \ 

puv 

pv 2 + p 

pvw 

V v(e + p) J 


and 


H = 


/ pw \ 

puw 

pvw 

pw 2 + p 
\w(e, + p) / 


In addition to eq. (2.7), the equation of state of an ideal gas is used: 


(2.9) 


p = pRT, (2.10) 

where p is the pressure, R the gas constant, and T the temperature. Since we assume that 
the incident wave hits the panel at t — 0, the flow on the transmitted side, y < 0, is at 
rest for t < 0. The initial data at t = 0 for y < 0 are 

u = v = w = 0 , p = Poo and p = Poo ■ (2-11) 

On the rigid plate we have zero normal velocity, 

u(a:,0,z) = 0 for {x,z) . (2-12) 

The motion of the flexible panel, V, is described by a system of three, nonlinear partial 
differential equations given by [10]; 

DV A y + pphpu + 7 y t = A p 

+ 2 [(«X + \nl){'Ux + VVzz) 

1 — v l 2 

+ ( w °z + \vl)( 7 hz + I'Pxx) 

+ (1 - i')v*z (a® + t a” + VxPz)} (2.13) 
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u° xx + diu° zz + d 2 w° xz = -7/ I (j/ XI + d X 7) zz ) - d 2 i]zVxz 
w° zz + di w xx + d 2 u° xz = -i] z (i] zz + d ] 7] xx ) - d 2 i] x i) xz 


where 


Vi] 1]xxxx + 2 J] XX zz + dz 

, _ 1 ~ " J _ 1 + J/ r»_ 

«1 = ~ 7 M 2 — 1 D — 


El t 3 


(2.14) 

(2.15) 


2 7 2 ’ _ 12(1 — // 2 ) ’ 
w° and u>° are the in-plane displacements, and 7/ is the transverse displacement. The 
physical constants for the panel appearing in (2.13) are; the stiffness ( D ), the density (p p ), 
the thickness (/i), the physical damping ( 7 ), the modulus of elasticity ( E ) and the Poisson 
ratio (//). The system of equations (2.13) is solved subject to the homogeneous initial 
condition 

at t — 0 u° = w° = r; = r/ t = 0 (2.1G) 


and the clamped boundary conditions on C given by 

X — Xq , X 0 + L , U° = W° — 7/ = 7] x = 0 , 

Z — Z 0 , Zo + W , U° = U7° = 7/ = 7/s = 0 . 


(2.17) 


In eq. (2.13), the load A p defined by (2.6) contains the coupling with the acoustic field p‘ 
and the forcing term (p'^ + ]/ ) which represents the load on the panel if the panel were 
rigid. We assume the forcing term to be a harmonic wave of the form 


p (0 = f .sin(wf) H(t ) on y = 0 + . 


(2.18) 


where e and u are the amplitude and frequency of the wave and H denotes the Heaviside 
unit step function. The load on the panel V given by (2.6) becomes 

A p = [ 2 p l — e .sm(u>t)] H{t) . (2.19) 


Another condition coupling the flow-fields and the panel oscillation is the kinematic con- 
ditions, 

7]t(t,x,z) = v ± (t,x,0 ± ,z) , ( x,z ) 6 V . (2.20) 
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They are imposed on y = 0"^ because of the small panel displacement, (1.1). The numerical 
scheme for the solution of the nonlinear plate equation (2.13) and the Euler equations (2.7) 
for the transmitted waves with homogeneous initial data will be described in Sec. 3. 

2.2 The Decoupled Model 

Now we analyse the nonlinear panel oscillation excited by a weak incident pressure 
wave, un der the assumption that the pressure fluctuation remains much smaller than 
the ambient pressure poo, i. e., 

(P - Poo)/Poo < 1 ■ (2.21) 

Consequently, the small disturbance theory is applicable to the flow fields in the incident 
and transmitted sides and the panel/acoustic interaction problem is described by (2.1) - 
(2.6). The velocity potential ^^(t, x, y, z) is governed by the simple wave equation, 

(C~ 2 dl - d\ x - d 2 yy - d] t ) ** = 0 , for ±y> 0 , (2.22) 

and the acoustic pressure and velocity v* are related to the potential by 

p ± = -p^dtS* and v ± = Vfc* . (2.23) 

In particular, we have v = d y $ where v denotes the vertical velocity. The vertical dis- 
placement of the panel, r;(t, .r, 2 ), is governed by the system of partial differential equations 
(2.13). To produce the special forcing term, (p (,) + p (r) ), acting on the panel specified by 
(2.18) in Sec. 2.1, the incident potential should be 

<t) (i) (t,x,y,z ) = — € — cos[<j(< + + ^) . (2.24) 

Note that the incident wave c.6 ^ 1 ■ is a solution of (2.22) in the whole space. When the panel 
is rigid, the incident wave is reflected by the rigid (x, z ) plane and the reflected wave is, 

<j) (r) {t,x,y,z) = <f>^ l \t,x,-y,z) for y> 0. (2.25) 
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On the side y < 0, the flow-field remains at rest because there is no transmitted wave. 

For a flexible panel, the panel oscillation excited by the pressure difference across the 
panel —2 p^(t, x, 0 + , z) induces a scattered wave <f>^ on the incident side and a transmitted 
wave (jy l \ Thus we write 


$ = ^(0 + tfr) + ^<*> , y > 0 (2.26a) 

and 

$ = jW , y < 0 . (2.266) 

Note that <J > ( s ■* and <f>^ are governed by the simple wave equation (2.22) in y > 0 and y < 0 
respectively. The kinematic conditions on the (x, z) plane, (2.20) and (2.23), become 

dy<j>( 3 \t, x,0 + ,z) = d t i](t,x,z) and d y <j) (l \t, x,0 _ , z) = d t 7j(t, x, z) , (2.27) 

where i / denotes the extension of the transverse displacement of the panel, V, to the (x, 
z) plane, i. e., 

r/ = 0 , (x,z) £ V . (2.28) 

If the incident wave front hits the panel at t = 0, we can impose the homogeneous initial 
conditions on (f >^ , <j>^ and r/, 

(j)^ = 0 , dt<f>^ = 0 , <j>^ = 0 , = 0 and y = 0 , d t y = 0 for t < 0 . 

(2.29) 

Since the ambient fluid above and below the panel are the same, (2.27) implies that the 
velocity potential induced by the panel oscillation has to be anti-symmetric in y, i. e., 

<j) {i) (t, x, y, z) = - , x, -y, z) , y < 0 . (2.30) 

Then the pressure difference across can be written as 

Ap = 2 Poo [d t <j> {i \t,x,0 + ,*) + d t <f> i3) (t,x,0 + ,z)] , (x,z) € V . (2.31) 
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The derivation of the system of equations, (2.21) to (2.31), completes the formulation of 
the structural/acoustic interaction problem using linear theory for the acoustic fields. The 
panel oscillation, which can be nonlinear is governed by the system of equations, (2.13) - 
(2.17). 

The velocity potential induced by the panel oscillation ?/ is governed by the wave 
equation (2.22), the initial conditions (2.29) and the boundary condition (2.27). The 
solution is given by the Kirchhoff formula, 


> 0,z) = J J 


R 


dx'dz' , 


(2.32) 


where R = [(x — x') 2 + y 2 + (z — 2') 2 ] 1 / 2 denotes the distance from a point (x, y > 0, z) to 
a source at (x 7 , 0, z 1 ) and {•} denotes the retarded value of, i. e., 


{r/ t (f,x',0,2 7 )} = T) t (t - ^,x',0,z') . 


(2.33) 


The domain of dependence of (j)( a \t,x,y, z) is the circular disc "H in the (x 7 , z') plane, i. 
e., 

H | R<Ct or r 2 = (x' - xf + (z 1 - z) 2 < C 2 t 2 - y 2 . (2.34) 

The domain of integration in (2.32) is the intersection of 7i and the panel, i. e., 


g = n n v . 


(2.35) 


Now we introduce the polar coordinates, r, 8 , centered at (x, 2 ), i. e., 


1 — x = r cos 8 and z' — z — r sin<? 


(2.36) 


and (2.32) becomes, 


= r,») ] ^ . (2 37) 


with 


g(t, r, 8) = r/(t, x + r cos 8, z + r sin 8) 
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(2.38) 



Here R = (r 2 + y 2 ) 1/f2 represents the slant height of a circular cone with vertex P(.r, y, z), 
a vertical axis and a base circle of radius r in the (a:, z ) plane. On account of (2.28), we 
can extend the domain of integration to 7i and rewrite (2.37) as an iterated integral, 

1 r('t p2n E? j* 

<f> (a) (t,x,y > 0,2) = [J de g t {t - g,r,Q) ] , (2.39) 

The integral in 8 represents the contribution of the sources on a base circle of radius r. 

In (2.31), we need to relate the unknown, <j>\ a \x, 0 + , z) on the panel to r/(t,.r,y) for 
(x, z) € V. This is obtained by differentiating the above equation with respect to t and 
using (2.28). This results in replacing ^ 3 ) on the left-hand-side of (2.39) and g% in the 
integrand by <// ( ^ and gtt, respectively. 

As y — » 0 + , R — * r and the circular disc 'H is bounded by the sonic circle v = Ct. 
Equation (2.39) yields 

, f2n rCt 

4 , \t,x,0 + ,z) = -±J I g u (t - dr M . (2.40) 

Note that in this form we remove the kernel 1 / R. in (2.32), which becomes singular as y — > 0 
and r — » 0 and shows that the area of the domain of integration is bounded above by 2irCt. 
By using (2.31) and (2.40), (2.13) becomes a system of integro-differrential equations for 
the panel oscillation, 7/. The initial and boundary conditions for t/ are (2.16) and (2.17). 
Thus we complete the formulation of the decoupled model. The numerical solution of this 
system will be described in Sec. 3. 

2.3 Application to Noise Transmission in a Supersonic Flight at High Altitude 

The analysis reported in Sec. 2.2 is based on two assumptions 

(i) The flow-fields above and below the panel have the same ambient condition, i. e., 
pto = PZo, pto = Pi o al l (1 = 0. 

(ii) The panel is planar and not prestressed, i. e., under zero tension. 
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The analysis can be employed to estimate the noise transmitted into the fuselage flying 
at low speed and a relatively low level, where the static pressure remains nearly equal to 
one atmosphere. Assumption (i) is valid when terms of the order of the flight Mach number 
M are omitted. Assumption (ii) is valid since the size of the panel L is much smaller than 
the radius of the fuselage, Rf. There is no hoop stress when under assumption 

(i). In practice, there is a pressure difference (p^ — p+>), even for a low Mach number 
flight. 

For an airplane cruising at subsonic or supersonic speed, it is necessary to account for 
the differences in densities and in pressures inside and outside the fuselage and the moving 
media outside and to include the contribution of the hoop stress to the panel oscillation. 
The formulation of the decoupled models for structural acoustic interactions in a moving 
media at subsonic, and supersonic speeds is carried out using the linear theory for the flow- 
field. The system of integro-differential equations for the panel oscillation is then reduced 
to a system of differential equations when the acoustic wave length is much smaller than 
the surface wave length [11]. 

In this subsection, we show that with minor modifications the above analysis becomes 
applicable to a supersonic airplane cruising at high altitude. We make use of the fact that 
the cabin pressure, p^, and temperature and hence the density are maintained near the 
ambient level while the outside pressure p+, and density p+, at the flight altitude H are 
much lower. The density ratio, p^/p^ is about 1/10 or 1/20 when H = 60,000 or 80,000 
ft. The same is true for the pressure ratio while the absolute temperature ratio and hence 
the ratio of the speed of sound remains 0(1), i. e., 

I‘ = pVpZ< 1, Pi/Pi = OM and C+/C- = 0( 1) . (2.41) 

Thus, the density ratio // serves as a small parameter in the following modification of the 
analysis in Sec. 2.2. 

On the transmitted side, y < 0, the governing equations, (2.13) (2.23), (2.26b), (2.27) 
and (2.29), remain valid when the superscript ” is added to poo and C. On the incident 
side, y > 0, equations (2.25), (2.26a), and (2.29) remain the same while (2.22), (2.23) 
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and (2.27) have to be modified to include the convection terms. Consequently, (2.30) and 
(2.31) will be modified. In particular, we note the acoustic pressure relationship (2.23) in 
the transmitted side, 

p~ -pZo =p (<) = , y< 0 (2.42) 

and that on the incident side, 

P + - pto = P (,) + P (r) + P {9) = -pp^[d t + MC + d x }$ , y > 0 . (2.43) 

Using (2.26a) and (2.43), we have 

P {s) = ~PP^[dt + MC + d x ]<f> {a) y > 0 , (2.44) 

Likewise, the kinematic, condition on the y = 0 + is modified and (2.27) becomes, 
d y <() i ' s \t,x,0 + ,z) — [d t + MC + d x ]r)(t,x,z) and dj,</>^(t, x, 0“, z) = d t r}(t, x,z ) . (2.45) 

With the convection term at most of the order of the unsteady term, or one order smaller 
under the short acoustic wave approximation, we deduce from (2.45) 

d y <f> w (t y x,0~,z) = O{d y <j> ia \t,x,0 + y z)) (2.46) 


and then from (2.44) 

p {a \t,x,0 + ,z) = 0(p p {t \t,x,0~,z)) . (2.47) 


By using (2.41) and (2.44), the load on the panel, which is the pressure difference, becomes 


A P = (pZ ~ rfo) + [-2p (0 + (P M - P M )W) 

= P~ [1 - 0(y)\ + [—2 p<‘> + p<‘>( 1 + 0(p)\H(t) . 


(2.48) 


The constant load (p^ — p+,) ~ p^ is present for all t. The unsteady load due to the 
acoustic waves is present only for t > 0 and is much smaller then p^. The constant load 
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produces a static deflection 7/o(^, ~ ) which does not generate far-field sound. The unsteady 
load excites panel oscillation ii'(t,x,z ) and generates the far-field sound. Therefore, we 
separate the static panel deflection from the unsteady oscillation, 

7/ = 7/o(x, z) + 7 /'(*, X, z ) . (2.49) 

Note that the panel oscillation, 7 /, is coupled with only the acoustic field (2.48) in the 
transmitted side, since p M is prescribed and p ^ is 0(p) related to This is equivalent 
to say that in a high altitude flight the acoustic, damping in the flow-field outside of the. 
fuselage is of the order of the density ratio p = ptJ PZo relative to that inside and hence 
is neglegible. 

Since the fluid is at rest inside the fuselage, the transmitted wave obeys the simple 
wave equation (2.22). We can make use of the analysis in Sec. 2.2, in particular, (2.30) 
and (2.40), to relate the on surface transmitted pressure to the panel oscillation, 

p w (t, x,0~,z) = — p^dt^ ~ — 7 ^ j j 9 dr d6 . (2.50) 

Thus we have a closed system for the panel oscillation, 7 /, from which we can then compute 
the transmitted and scattered acoustic fields. 

3. Numerical Methods 

For the fully coupled model , the nonlinear Euler equations, (2.7), are solved using an 
explicit finite difference scheme. The scheme, which is a generalization of MacCormack s 
scheme obtained by Gottlieb and Turkel [12], is fourth-order accurate in space and second- 
order accurate in time. Further details on the implementation of the scheme can be found 
in Frendi et al. [1,2]. The physical boundary of the computational domain (the bottom 
boundary, see Fig. 1) is composed of a flexible panel clamped between rigid plates. Over 
the rigid plates, the vertical velocity is zero (v - 0) and the surface temperature is T w , 
which in this paper is the same as the fluid temperature, T 0 c . The x and 2 components 
of the velocity ( u and to, respectively) are obtained through linear extrapolation from the 



interior of the computational domain. The pressure boundary condition is obtained using 
the normal momentum equation by simply imposing the normal gradient of the sum of 
pressure and vertical momentum flux to be zero, i. e., 0 y (p pv 2 ) = 0. Over the flexible 
panel, the vertical velocity is set to be equal to that of the panel and the temperature is 
T w . The x and 2 velocity components and the pressure are extrapolated from the interior 
of the computational domain. The appropriate boundary conditions on the computational 
domain are derived using the method of characteristics [13]. One should mention that 
the characteristic boundary conditions and the extrapolation are only first order accurate, 
while the interior scheme is fourth order. These boundaries are believed to be a source for 
numerical error in the fully coupled model. The dimensions of the computational domain 
are, 61 cm, 3.05 m and 40.64 cm in the x, y and 2 directions, respectively. The number of 
computational points used are, 121, 241, and 81 in the respective directions. 

The nonlinear plate equations, (2.13), are solved using a finite element method devel- 
oped by Robinson [14]. The panel is 30.5 cm long, 20.32 cm wide and 0.102 cm thick, and 
the number of elements used are 6 and 8 respectively. 

Since the grid used for solving the plate equations is rectangular, it is easlier to 
evaluate the integral in (2.40) in cartesian coordinates. In this case there is a singularity 
at each computational point corresponding to r — 0. In order to overcome this difficulty, 
a Taylor series expansion up to the second order is used to calculate the contribution of 
the singular point. The contribution of the various points on the panel that lie within 
the sonic circle ( R < Ct) is calculated by first integrating in x using a combination of 
Simpson and trapezoidal rules. The result is then integrated in z using Simpson’s rule 
of integration. Because of the presence of the retarded time in the integral, the vertical 
velocity of the plate (?/<) is stored at each point for several time-steps. The evaluation 
of the double integral can use a much larger time step than the At for the integration of 
the differential equations, yet having the same degree of accuracy. For a plate of given 
dimensions (L, IF), and for a fixed time-step At, the maximum number of time-steps to 
be stored is N = \/L 2 -\- W 2 j (C At). The number N has to be changed when calculating 
the radiated pressure away from the plate. The number N is related to the radius of the 
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largest sonic circle in the domain that contains the flexible panel. For the cases presented 
in the next section N = 1503 for a time-step At = 4xl0 -6 sec. 

4. Results and Discussion 

The numerical schemes for the two models presented in Sec. 3 are used to predict the 
vibration of a flexible panel and the resulting acoustic radiation. The panel is forced to 
vibrate by harmonic plane acoustic waves at normal incidence. The frequency of the source 
is 751 Hz, which corresponds to a natural frequency of the panel. Two different amplitudes 
of the incident waves are used. The properties of the panel, which are considered to be 
uniform, are: density pp = 4450.15 Kg/m 3 , modulus of elasticity E = 1.10316x10 N/m , 
Poisson ratio u = 0.33 and a damping ratio of 0.01. The acoustic fluid properties are: 
temperature Too = 288.33 °K, density Poo = 1.23 Kg/m 3 , pressure p 0 Q = 1.013x10 N/m 
and sound speed = 340 m/sec. The specific heat at constant volume is c v = 1.004 
KJ /(Kg °K), the ratio of specific heats is 7 = c p /c v = 1.4. In the-far field, the fluid is at 
rest. The variables plotted on figures 2-5 are nondimensional. The reference quantities are 
given by 


(X, y, Z , 1]) re f = Iref, *re/ = , 

drj 


C 2 

and T re f — — 


pre f = Poo, (u,V,W,-^) re f = C, and (P, c) re / - pooC . 

where the reference length is l re f = 0.3048 in. 


(4.1) 


Figure 1 shows the configuration of the computational domain, a rectangular box 
with the lower side composed of a flexible panel clamped between rigid plates. For a low 
excitation amplitude, 100 dB or 5 x 10 ^ atm., Fig. 2 shows that the panel response is 
linear. In this case, both models predict the same panel response as shown by the figure. 
This result is expected since the incident wave is extremely weak. Figure 3a-b show the 
time histories of the radiated near- and far-field pressures, 2.54 cm and 1.524 in, or L/ 12 
and 5 L, away from the panel center respectively. The radiated pressure predicted by the 
linear theory (2.32) is in excellent agreement with that predicted by the Euler equations 

(2.7). 
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When the level of the excitation is increased to 160 dB, or 0.01 atm., the response 
of the panel becomes nonlinear as shown by Fig. 4. The time history of the panel center 
displacement shows a non-periodic behaviour characteristic of nonlinear response. The 
figure also shows that even when the panel oscillation is nonlinear, the predictions of the 
two models are in reasonably good agreement, because the pressure variation is still much 
smaller than poo = 1 atm. The radiated near- and far-field pressure time histories are 
shown on Fig. 5a-b. In the near-field, Lf 12, the predictions of the two models are in 
reasonably good agreement and the maximum pressure variation is of the order of 0.02 
atm. In the far-field, 5 L, as shown in Fig. 5b, the maximum pressure variation is of the 
order of 10~ 4 atm. and is reduced by a factor of 100 from that in the near-field. In Fig. 
5b, the difference between the “decoupled model” and the “fully coupled model” can be 
obseved. But this difference is of the order of 1 0 times the pressure variation in the 
near-field and is within the accuracy of the numerical solution. 

From a computational view point, it is important to compare the performance of the 
two models based on the CPU time required by each calculation. In the linear vibration 
regime, the “fully coupled model” used 36000 seconds of CPU time on a Cray-ymp to 
advance the calculation by 10000 time-steps, whereas the “decoupled model” used only 
1000 seconds for the same calculation. In the nonlinear vibration regime, grid refinements 
were needed to resolve the large gradients both on the panel and in the radiation field. 
Therefore, in order to advance the calculation by 10000 time-steps, 72000 seconds were 
used by the “fully coupled model” while the “decoupled model” used only 1200 seconds. 

5. Conclusions 

An efficient model for coupling the vibration of a panel to the on surface acoustic radi- 
ation was derived. The model uncouples the panel vibrations from the acoustic wave prop- 
agation problem. The results showed that this model, referred to as “decoupled model”, 
accurately predicts the panel response and acoustic radiation in the linear and nonlinear 
vibration regimes so long as the pressure variation in the flow-field remains much smaller 
than the ambient pressure. For the cases studied in this paper, the computational cost of 
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the numerical integration (if this model is 36 times cheaper in the linear regime and 60 
times cheaper in the nonlinear regime than the cost of the “fully coupled model”. 
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Figure 1: Computational Domain 
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Figure 2: Time history of the displacement of the panel center for an excitation 
amplitude of 100 dB (linear). 
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Figure 3: Time history of the radiated pressure along a line perpendicular to 
the panel center, (a) 1 inch away from the panel, (b) 60 inches away from the 
panel. 
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Figure 5: Time history of the radiated pressure along a line perpendicular to 
the panel center, (a) 1 inch away from the panel, (b) 60 inches away from the 
panel. *U.S. GOVERNMENT PRINTING OFFICE: 1993-728-150/60050 
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